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Abstract 

We formulate and prove a periodic analog of Maxwell’s theorem relating stressed planar frameworks and their 
liftings to polyhedral surfaces with spherical topology. We use our lifting theorem to prove deformation and rigidity- 
theoretic properties for planar periodic pseudo-triangulations, generalizing features known for their finite counter¬ 
parts. These properties are then applied to questions originating in mathematical crystallography and materials sci¬ 
ence, concerning planar periodic auxetic structures and ultrarigid periodic frameworks. 

Keywords: periodic framework Maxwell’s theorem periodic stress liftings periodic pseudo-triangulation expansive 
motion auxetics ultrarigidity 


1 Introduction 

A remarkable correspondence between planar stressed graphs, their duals and polyhedral surfaces with a spherical 
topology has been established in 1870 by James Clerk Maxwell: 

Maxwell’s Theorem [36] A planar geometric graph (G^p) supports a non¬ 
trivial stress on its edges if and only if it has a dual reciprocal diagram 
and, at the same time, if and only if it has a non-trivial lifting to 3D as a 
polyhedral terrain. 

The necessary definitions are recalled below in Section 2. A closely re¬ 
lated instance of this theorem is the classical duality between Voronoi dia¬ 
grams and Delaunay tesselations, where the 3D lifting is onto a paraboloid. 
Maxwell’s diagrams, further popularized in Cremona’s book [16], were 
widely used for engineering calculations throughout the 19th and 20th cen¬ 
turies. The theorem has many other applications, for example in problems 
of robustness for geometric algorithms, rigidity theory, polyhedral combi¬ 
natorics and computational geometry [15, 30, 43, 44, 12, 48]. Most relevant 
to our undertaking is its role in establishing the existence of planar expansive motions used in the solution to the Car¬ 
penter’s Rule problem [12], and in proving the expansive properties of pointed pseudo-triangulation mechanisms that 
are central to the algorithm for convexifying simple planar polygons of [47, 48]. 

Our results. In this paper we prove the following periodic analog of Maxwell’s theorem. Figure 2 illustrates the 
concepts, whose precise definitions will be given in Section 3. 

Main Theorem Let (G,r,p, n) be a planar non-crossing periodic framework. A stress induced by a periodic lifting 
is a periodic stress and conversely, any periodic stress is induced by a periodic lifting, determined up to an arbitrary 
additive constant. 

Non-crossing periodic graphs can be seen as graphs embedded on the flat torus. However, as it will become clear 
from this paper, to reason on a fixed torus would be too restrictive a perspective. The most important foundational 
element that makes the new result possible is our recent deformation theory of periodic frameworks [3], which allows 
the periodicity lattice to deform. The corresponding notion of periodic stress is precisely the notion of stress that 
is needed for the Main Theorem. This stress is more constrained than the direct generalization of the classical self¬ 
stress used for finite frameworks, which is based solely on equilibrium at ah vertices. In order to maintain the proper 
distinction, we refer to the latter type of stress as an equilibrium stress. 



Figure 1: A finite planar stressed graph 
and a Maxwell lifting. 
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Figure 2: The connection between periodic stresses and liftings proven in this paper, (a) A stressed periodic framework. 
Vertex and edge orbits are identically colored, and a fundamental polygon of the periodicity lattice is highlighted in 
grey, (b) Coloring the faces helps visualize the 3D lifting of the framework as a periodic arrangement of “cubes”, (c) 
The same periodic graph, in a placement with no non-trivial periodic stresses, (d) The face coloring visually confirms 
that this is not the projection of a polyhedral surface: the faces do not look “fiat” in 3D. 


Our theorem was motivated by questions arising in mathematical crystallography and computational materials science. 
We demonstrate its usefulness with two applications: constructions of ultrarigid frameworks and auxetic mechanisms. 

Ultrarigidity of periodic frameworks. Leaving some 
technical details aside, our proof of the correspondence 
between periodic liftings and periodic stresses will pro¬ 
ceed by showing how to obtain a transparent, algebraic 
matching of all the concepts involved after a sufficient 
relaxation of periodicity . Relaxations of periodicity are 
a central concern in displacive phase transitions [20, 6] 
and successive relaxations give rise to difficult and im¬ 
portant problems for estimating the asymptotic behavior 
of a periodic framework. 

By definition, a periodic framework is ultrarigid if it is 
and remains infinitesimally periodically rigid under arbitrary relaxations of periodicity to subgroups of finite index (see 
Figure 3). This concept was introduced in [5], and illustrated with a few examples of crystalline materials exhibiting 
this property. Ultrarigidity provides a rigorous tool for studying the asymptotic rigidity of a periodic framework 
when successive relaxations of periodicity are applied. The proof techniques developed in this paper will allow us to 
construct infinite families of ultrarigid examples. 

For some related considerations and observations on rigidity and relaxation, we 
mention [41, 13]. The basic theory of periodic frameworks from the point of view of 
rigidity and fiexibility can be found in [3, 4, 6]. For wider or complementary aspects 
of periodic framework theory we suggest [50, 11] and references therein. 

Periodic pseudo-triangulations. We use the Main Theorem to study a new 
class of planar non-crossing periodic frameworks called periodic pointed pseudo¬ 
triangulations or shortly periodic pseudo-triangulations. They represent a natural 
analog of the finite pointed pseudo-triangulation frameworks defined and studied in 
[47, 48] and possess mutatis mutandis many outstanding characteristics related to 
rigidity and deformations [47,48,45,46]. Here we focus on the expansive one-degree- 
of-freedom mechanisms they provide and on the property of turning into ultrarigid 
periodic frameworks after the insertion of a single (adequate) edge-orbit. 

Deformations of periodic frameworks: auxetic and expansive behavior. The sig¬ 
nificance of the idea of expansive motion is well recognized in the finite setting [12, 48,46]: when the distance between 
any pair of vertices cannot decrease, self-collision of the framework is avoided. In the periodic setting, expansive 



Figure 4: A periodic pointed 
pseudo-triangulation. 



Figure 3: A periodic framework and a 2 x 1 relaxation of 
its periodicity lattice. 
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mechanisms have not been explicitly considered before, although a related, yet weaker notion of auxetic behavior has 
recently attracted a lot of attention in materials science [21, 33, 39]. Since the existing literature on auxetics is based 
on elasticity theory, we include a brief and necessarily selective overview of the relations existing between the purely 
geometric theory pursued in this paper and the larger context of periodic structures explored in crystallography, solid 

Interest in crystal morphology and structure motivated mathemati¬ 
cal studies of symmetry, lattice sphere packings and crystallographic 
groups [14, 34]. With the advent of X-ray diffraction, materials 
science gained access to atomic-scale configurations and bonding 
networks. An explicit mention of a periodic framework deforma¬ 
tion appears in Pauling’s 1930 paper [40]. Nevertheless, such ge¬ 
ometric investigations addressed only specific crystalline materials 
and remained mostly concerned with a number of instances of one- 
parameter deformations related to particular displacive phase transi¬ 
tions, as in [18, 20]. 

The notion of auxetic behavior is formulated using the concept of 
negative Poisson’s ratio [22, 23, 21], which relies on physical prop¬ 
erties of the material: when two forces pull in opposite directions 
along an axis, most materials are expected to expand along this axis and to contract along directions perpendicular 
to it. Auxetic behavior refers to the rather counter-intuitive lateral widening upon application of a longitudinal ten¬ 
sile strain. A purely geometric expression of this behavior is not anticipated in all situations. However, for periodic 
frameworks, we have recently proposed the general geometric notion of auxetic path in the deformation space of the 
periodic framework [8]. Relying on this formulation, we prove that an expansive deformation path is necessarily 
an auxetic path. Periodic pseudo-triangulations thus exhibit auxetic behavior and offer an infinite supply of planar 
examples of “auxetic frameworks”. By contrast, only a limited collection of sporadic and artisanal auxetic periodic 
examples (Figure 5) has appeared in the literature. One may already notice that the reentrant framework from Figure 5 
looks “almost” like a pointed pseudo-triangulation, except that it is not maximal. It needs one more edge-orbit for all 
its faces to become pseudo-triangles, as illustrated in Figure 16 of Section 6. 

Organization of the paper. In Section 2 we define the basic concepts needed to study the correspondence between 
liftings and stresses in (finite or infinite) frameworks. Section 3 specializes these concepts to periodic liftings and 
stresses. Section 4 concludes the proof of the Main Theorem by providing the necessary link with periodic rigid¬ 
ity and fiexibility. The Main theorem is applied in Section 5 to prove the expansive properties of periodic pseudo¬ 
triangulations. The connection with auxetic behavior and the ultrarigid character of periodic frameworks obtained 
from pseudo-triangulations are presented in the final Sections 6 and 7. 

2 Liftings and stresses 

To formulate and prove our Main Theorem, we start with those concepts and properties that do not depend (yet) on 
periodicity, which is introduced in the next section. 

Graphs and frameworks. We consider finite or countably infinite graphs which are simple (i.e. without loops or 
multiple edges), unoriented and of finite degree (or valency) at each vertex. Such a graph is given as a pair G = (V, E), 
with V the set of vertices and E the set of (unoriented) edges. We use lower case symbols w, v, • • • for vertices in V. 
An edge e = {w,v} = {v,w} G E has two endpoints u^v ^\/ ^uf^v and can be given two orientations (w,v) and (v,w). 
Two edges , ^2 ^ E are adjacent if they have a common endpoint: |^i fl ^ 2 ! = 1- The set of edges incident to a vertex 
V consists of all edges ^ G E having v as one of their endpoints. The degree (or valency) of a vertex is the number of 
edges incident to it. 

A placement of G in is given by a mapping p : \/ \-^ of the vertices to points in such that the two endpoints 
of each edge e = {u^v} e E are mapped to distinct points in W^\ p{u) 7 ^ p{v). An edge {u^v} is seen geometrically as 
an edge-segment [p{u)^p{v)], and an oriented edge {u^v) determines an edge-vector p(y) — p{u) G W^. We work under 


state physics and materials science [34, 50, 22]. 



Figure 5: This so-called ‘reentrant’ structure of 
hexagons is often used to illustrate auxetic be¬ 
havior. 
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the assumption that all placements are locally finite maps, that is, the preimage of any bounded set is finite. This is 
certainly true for the periodic placements defined in the next section. 

A framework or geometric graph (G^p) is a graph G together with a placement p, restricted in this paper to M? or M?. 
We use the term planar placement for M?, when the distinction is necessary. 


Planar non-crossing frameworks. A planar placement is non-crossing if any pair of edges induces disjoint closed 
segments, with the possible exception of the common endpoint, in the case when the edges are adjacent. A graph G 
is planar^ if it admits a non-crossing placement {G,p). We consider only connected graphs, therefore a non-crossing 
placement (G,/?) induces a connected subset of the plane, made of the points p(v),v G V and the edge segments 
[p{u),p{v)] with {m,v} G E. a face G is a connected component of the complement of the placement (G,p). We 
assume throughout the paper that the boundary of each face is a simple finite polygon. A face is described 
combinatorially by the cyclic collection of its boundary vertices or edges, and each edge is on the boundary of exactly 
two distinct faces. 

For a planar non-crossing framework it is convenient to use now the same symbol G to denote the entire collection 
G = (V, E, F) of vertices V, edges E and faces F. For notational simplicity, we allow the capital letters G,y,... to stand 
for open face domains, closed faces, or the corresponding boundary cycles (as needed in various contexts). 

When referring to a planar graph G (with no reference to a particular placement) we assume that the choice of face 
cycles F is also given, i.e. G = (V, E, F). We note that even when the underlying graph G = (V^E^F) of a framework 
(G^p) is a planar graph, the particular placement p of the framework may have crossings: we still refer to the 
realization of a face, although it may be a self-interseating polygon. 

To a planar graph G = (V, E, F) we can associate a dual structure G* = (V*, E*, F*) defined as the abstract triple whose 
vertices V* correspond to the faces F of G i.e. V* = F, and whose edges E* are in one-to-one correspondence with the 
edges E of G, as follows: if two faces U and W share an edge e, then the dual vertices U* and W* are connected by 
the dual edge e*. The dual faces F* correspond to the vertices of G i.e. F* = V, with the cycle of faces around a vertex 
inducing its corresponding dual face. By abuse of language, we may refer to G* = (V*, E*, F*) as the dual graph of 
G = (V, E, F), although it may have multiple edges. 


Orientation rule. In a planar non-crossing framework, an edge 
{m, v} induces a segment [p{u),p{v)], and it belongs to the boundary 
of exactly two faces, say U and V. In the dual graph G* = (F, E) 
these two faces U and V represent two vertices connected by the un¬ 
oriented edge {G, V } dual to {w, v}. Later on, we will need to match 
an oriented edge in the primal graph G = (V, E) with an orientation 
of its dual edge in G* = (F, E). We use the following convention. 
The oriented edge segment [p{u)^p{v)] gives opposite senses for go¬ 
ing around face G and face V, say counter-clockwise around V and 
clockwise around G. Then the matching orientation of the edge in 
the dual graph G* is from G to V. In short: ‘from clockwise to 
counter-clockwise’, as illustrated in Figure 6. 

The matching orientation described above for elements in the edge 
set E, when considered as oriented edges in G = (V, E), respectively 
Figure 6: The orientation convention for face G* = (F, E), gives a well-formed double pair {{u,v), (G,y)), or sim- 
edges and vertex edges. ply a tetrad (w, v, G, y ). In computational geometry, these tetrads are 

implicit in the quad-edge data structure used for representing general 
surfaces [29]. Reversing orientation on the edge gives the tetrad (v, u^V^U). Below, we will refer to cycles (cdllod face- 
cycles) of oriented edges in the dual graph G* = {F^E): the orientation rule described above gives an unambiguous 
correspondence with oriented edges {u^v) in the primal graph G and their corresponding edge vectors p{v) — p{u) in 
a geometric placement p of G. 

^Note that we use planar for the graph, as is customary in graph theory, and non-crossing for the framework. Our use of planar framework is 
customary in rigidity theory, and refers to a placement in the plane. 
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Stressed frameworks. An equilibrium stress or, shortly, a stress on a planar (finite or infinite, possibly crossing) 
framework^ is an assignment s \ E of scalar values {se}eeE to the edges £" of G in such a way that the edge 
vectors incident to each vertex u eV, scaled by their corresponding stresses, are in equilibrium, i.e. sum up to zero: 

^ Se{p{v) — p{u)) = 0, for fixed j/gV (1) 

e={u,v}EE 

When all Se are zero, the stress 5 ' is trivial, when all Se 7 ^ 0, the stress is called nowhere zero. The space of all 
equilibrium stresses of a framework is a vector space, so if a framework has a non-trivial stress, then it is not unique; 
in particular, any rescaling of it is also a stress. 

Lifting. A lifting of the planar framework (G,/?) is a continuous function H whose restriction to any face 

is an affine function. The lifting assigns a height, or altitude H{q) to each point q in R^ (seen as the plane z = 0 in 
R^), in such a way that the lifted faces are fiat (all face cycle vertices lie in the same plane) and connect continuously 
along the lifted edge segments. The height function is completely determined by the values H{p(y)) at the vertices of 
the framework, and its graph appears as a polyhedral surface or terrain over the face-tiling in the reference plane. A 
lifting is trivial if all its faces are lifted in the same plane, that is, when H is affine on A lifting is strict if no two 
adjacent faces are lifted to the same plane. 

With these concepts in place, we move on to the correspondences involved in Maxwell’s theorem. 

Equilibrium stress associated to a lifting. Let ^ be a lifting of a framework (G, p). With usual dot product notation, 
the expression of H restricted to a face U takes the form 

H{q) = Vu-q + Cu, fovqeUd^ (2) 

where Vu € is the projection on the reference plane of the normal to face U and Cy G M. 

The system of vectors and constants H = {Vu^Cu)uef is subject to the compatibility conditions on edges {u^v} shared 
by pairs of adjacent faces {U^V}: 


Vu • p{u) -\-Cu = Vv ' p{u) + Cy 

V[/ • p(y) + C[/ = Vy • p{v) + Cy (3) 


From this, we infer that the vector Vy — Vu is orthogonal to the edge vector p{v) — p{u). Fig. 7 illustrates the relation¬ 
ship. 



Figure 7: The normals to two adjacent 
lifted faces induce the dual orthogonal 
edge. 


Given a tetrad (t/,v,f/,y) of dual edges, and using the notation (v,y)^ = 
(—y, v) for the clockwise rotation with ;r /2 of a vector (v,y) G we define 

the stress factor on the edge {u^ v} G E, associated to the lifting H, as being 
the proportionality factor Suv given by: 


W-Vu = Suv{p{v) - p ( m ))^ (4) 

Since the sum involves the vectors around a closed polygon (the face-cycle 
around a vertex), the equilibrium condition (1) is satisfied. This proves: 

Proposition 1. Eor any lifting H of a planar non-crossing framework 
(G^p), there exists a canonically associated equilibrium stress on the 
framework. 

The correspondence between liftings and stresses described above is essen¬ 
tially the one given by Maxwell, who formulated it through the following 
geometric construction. The normal direction to the planar region corre¬ 
sponding to a face G G F in the lifted terrain is given hy Nu = (V[/, — 1) G 


^Also called a self-stress in the rigidity theory literature. 
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U ef. When all these normal vectors are taken through the point (0,0,1), they intersect the reference plane z = 0 
in the system of points {V[/ G The classical “theorem of the three perpendiculars” implies the orthogonality 

(vy — Vu) • {p{v) — p{u)) = 0 observed above. In Figure 7 we see the normals to two adjacent lifted faces, taken from 
a point of the lifted common edge. They intersect the reference plane in N{U) and NiV), with vectors p{v) — p{u) and 
N{V) -N{U) orthogonal. 

Reciprocal diagram. A framework (G*,p*) associated to the dual graph G* of a planar framework (G,/?) is called 
a reciprocal diagram if the corresponding primal-dual edges are perpendicular. If in the previous construction we 
join the points {V[/ G | U G F} by edges dual to the primal ones, we obtain a reciprocal diagram associated to 
the lifting H. We note that it is possible for several vertices Vu to coincide, and this happens precisely when several 
planar regions in the lifting have identical normal directions. An extreme case arises for liftings with globally affine 
functions H. They give a planar (trivial) terrain over the reference plane, have constant Vu and induce the trivial stress 

{Se}e^E = 0 . 

From equilibrium stresses to liftings. The direction from stresses to Maxwell liftings requires more work. 

Proposition 2. Let s = (*s’g)eeE be an equilibrium stress for the framework (G^p). Then there exists a lifting H which 
induces s, determined up to addition of a global affine function. 


Proof We have to find a set of parameters (Vu^Cu), indexed by faces and satisfying the conditions (3) and (4) in terms 
of the given placement p. Let us choose an initial face Go with an arbitrary lifting {Vuq^Cuq) = (Vo,Co). We show that 
once this initial choice has been made, the lifting is then uniquely determined by the placement p and the stress values 

5 '. 


For this purpose, we solve the linear system (3) in a step-by-step manner, progressing from face to adjacent face, 
starting at Gq. We remark that an infinite graph induces an infinite such linear system, but nevertheless it can be solved 
incrementally, as we now show. We consider a path through adjacent faces labeled Go, Gi,..., G„, with corresponding 
liftings {Vi^Ci) and successive tetrads (/7/,^/,G;,G/+i). Thus, [pi^qt] is the common edge between faces Ui and G^+i, 
with the proper orientation induced by the direction in which we walk through the path of adjacent faces, and the 
orientation rule described earlier. The given stress on this edge is denoted here by Si. 

We first compute the parameters Vk by unfolding the relationships (4) along the path from Go to Uk'. 

k 

Vk+i = yk+Sk{qk-Pk)^ =... = + ( 5 ) 

i=0 

Similarly, we compute the parameters by unfolding the first of the two relationships (3) along the path from Go to 
Uk, and using (4) at each step: 


k 

Q+1 = Q - (Vy^+I -^k)' Pk = = Gq - 'Y^sfqi- Pi)^ 'Pi (6) 

i=0 

Using the identity (qf—pi)^ • Pi = det{qi pt) and the notation \qi pi\ := det{qi pt) in (5) and (6), the expression (2) of 
the height function becomes: 

n—\ n—\ 

H{q) = Vn-q + Cn = {Vo+Y,Si{qi-pi)^)-q+{Co- Y, Pi\) ( 7 ) 

i=0 i=0 

It remains to check that the expression (7) is independent of the face-path chosen from Go to G„. To verify this property, 
we have to check that the following sums vanish for any face-cycle: 

Y Si{qi-pi)=0 and Y Si\qiPi\=Q (8) 

i^f ace—cycle i^f ace—cycle 
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Figure 8: (a) A Delaunay framework (in black) corresponding to a periodic set with four site orbits under maximal 
translational symmetry. Dual vertices are shown as colored centers of the faces. Primal-dual edge pairs are orthogonal. 

(b) The equilibrium stress induced by the reciprocal diagram has a non-periodic lifting to the paraboloid z = 

(c) The reciprocal diagram of the periodic graph on the left is the Voronoi diagram of the periodic sites. 


It suffices to verify these relations over face-cycles corresponding to simple topological loops. In this case, Jordan’s 
simple curve theorem gives a set of vertices inside the loop. When we sum the equilibrium stress condition (1) over 
the vertices inside the loop, terms cancel in pairs for adjacent vertices and yield the first identity. For the second 
part, we observe that the equilibrium stress condition (1) induces the identity |L{M,v}eE‘^Mv(p(v) — p{u)) p{u)\ = 
0, for fixed vertex u. Rewritting it as L{M,v}eE ^uv\p{^) p{u) \ = 0 and again taking the sum over the vertices inside the 
loop, we obtain the second identity in (8). 

Since the initial choice (vo,Co) was arbitrary, the lifting H is determined only up to a global affine function. □ 

Mountain-valley edge liftings and stress factor signs. We review now an elementary geometrical criterion for 
determining the sign of a stress factor Suv in a stress 5' associated to a lifting. Let us consider the lifted 3D polygonal 
surface projecting over the faces of the graph in the reference plane, as in Figure 1, and let uv be part of a tetrad 
(w, V, A vertical translation of the whole polygonal surface does not affect the stress considerations. We assume 
therefore that the lifting of the two faces U and V is above the reference plane. Figure 7 illustrates the case of a 
negative stress factor Suv, corresponding visually with a lifted edge that looks like a ‘mountain ridge’. The normal 
directions to the lifted faces over U and V are here taken through a point of the lifted common edge and intersect the 
reference plane in points denoted N{U) and N{y). For the tetrad (t/, v, f/, V), the sign of the stress factor is that of the 
determinant |;?(v) — p{u) N{V) —N{U) \ and is consistent with the orientation rule described in Section 2. 

For the geometrical interpretation, let us look at the dihedral angle of the lifted faces, understood as the intersection of 
the two half-spaces below the respective supporting planes. Then, when the measure of this dihedral angle is less than 
7t, we have a negative stress factor and the visual landscape suggestive of a mountain range, while for a dihedral angle 
greater than n, we have a positive stress factor and the visual landscape suggestive of a valley. When the dihedral 
angle is exactly n, the stress factor is zero and the two lifted faces are in the same plane. A non-fiat edge in the lifting 
is therefore said to be a mountain edge if the terrain is concave in its neighborhood, and a valley otherwise. The 
correspondence between stresses and liftings is now refined by this well-known property [15]: 

Proposition 3. The Maxwell correspondence between stressed graphs and liftings takes planar edges with a negative 
stress to mountain edges in the 3D lifting, those with positive stress to valley edges and those with zero stress to the 
common plane of the two adjacent lifted faces. 

Comment. As mentioned, this correspondence between liftings and stresses is a direct adaptation of the one for¬ 
mulated by Maxwell in the finite setting [35, 36] (see also the expositions given in [30, 15, 43]). The arrangement 
and presentation given in this section, including the notation, are meant to serve as a preamble for the more elabo¬ 
rate periodic version, which relies on notions developed recently in our deformation theory of periodic frameworks 
[3,4, 6]. 

We conclude this section with a noteworthy example, illustrated in Fig. 8, which will be used in the next two sections 
to illustrate the critical distinction between equilibrium stress and periodic stress. 
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Example 1. (Delaunay tesselations and Voronoi diagrams of periodic point sets) The classical Voronoi-Delauney 
duality, applied to a (countable locally finite) periodic point set p in yields a dual pair of non-crossing periodic 

frameworks whose corresponding dual edges are orthogonal. For a fixed point pi in the given set p, its (open) Voronoi 
cell is defined as the planar domain made of all the points in M? whose distance to pi is strictly less than the distance 
to any other point in the set p. The union of all these cell boundaries give the vertices and edges of the Voronoi 
framework. The Delauney framework is its dual, with vertices at the original point set p and with edges orthogonal 
to those of the Voronoi framework; its faces are convex polygons whose vertices are cocircular, and the circumscribed 
dries are empty of any other vertices in the given point set [17]. Fig. 8 gives an illustration with periodicity lattice 
generated by the standard basis. There are four orbits of vertices, with representatives at the corners of one of the 
diamonds in Fig. 8(a). All faces are inscribed polygons and the centers of the corresponding circumscribed circles 
give the vertices of the Voronoi diagram. 


3 Equilibrium and periodic stresses on periodic frameworks 

We turn now to the main object studied in this paper, the infinite periodic framework as defined in [3, 4]. We emphasize 
from the outset that a periodic graph is not just an infinite graph, as it was in the previous section: the definition includes 
a periodicity group of graph automorphisms. In this paper we focus on the planar case d = 2 and on connected non¬ 
crossing periodic frameworks. Some statements and constructions will be valid in broader contexts, but for the sake 
of a streamlined presentation we will stay within this class. 

Planar periodic frameworks. A 2-periodic framework [3, 4], denoted as (G,r,/7, ;r), is given by an infinite graph 
G, a periodicity group F acting on G, a placement p and a representation n. The graph G = (V, E) is simple (has 
no multiple edges and no loops) and connected, with an infinite set of vertices V and (unoriented) edges E. The 
periodicity group F C Aut{G) is a free Abelian group of rank two acting on G without fixed vertices or fixed edges. 
We consider only the case where the quotient multigraph G/F (which may have loops and multiple edges) is finite, and 
use n= I V/F| and m = | E/F| to denote the number of vertex and edge orbits. The function p'.V gives a specific 
placement of the vertices as points in the plane, in such a way that any two vertices joined by an edge in E are mapped 
to distinct points. The injective group morphism 7t '.T ^ gives a faithful representation of F by a lattice of 

translations ;r(F) = A of rank two in the group of planar translations £7 (M^) = Mf. The placement is periodic in the 
obvious sense that the abstract action of the periodicity group F is replicated by the action of the periodicity lattice 
A = ;r(F) on the placed vertices: p{yv) = n{y){p{v))^ for all 7 G F, v G V. 

Non-crossing planar periodic frameworks G = ((V, E, F),F) in They have an underlying planar graph G = 
(V, E,F) with the natural action of the periodicity group F. The dual 2-periodic graph G* = ((V*, E*, F*),F) is 
obtained from the abstract dual of the infinite graph G = (V, E, F), with the periodicity group F acting on it in the same 
manner as it acts on the sets of the primal graph F = V*, E = E*, V = F*. If we denote by = card{V/T) the number 
of face orbits under F, then Euler’s formula for the torus M^/A D G/F gives the relation: 

^ —m + ^* = 0, that is n-\-n*=m. (9) 


F-invariant equilibrium stress on a periodic framework. An equilibrium stress 3’ of the periodic framework 
(G,/7,F, ;r) is called a F-invariant equilibrium stress if it is invariant on edge orbits E/F. A F-invariant equilib¬ 
rium stress can be calculated by solving a finite linear system of equations of type (1), where the unknowns are the 
stresses for the edge representatives in E/F and the equations correspond to equilibrium conditions for the vertex 
representatives in V/F. 

Preview: periodic stress from periodic lifting. Not all equilibrium stresses of a periodic framework are F-invariant. 
Although we have not yet defined periodic stresses (Definition (6) in the next Section 4), we mention, as a preview, 
that they must be F-invariant. But this will not be sufficient: for some periodic frameworks, the periodic stresses are 
a strict subset of the F-invariant ones. Figure 9 illustrates such a situation, using colors to indicate the stress sign (red 
for positive, blue for negative stress, gray for zero stress) for the various types of stresses that a periodic framework 
may support (not F-invariant, F-invariant non-periodic and periodic). 
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Figure 9: The Voronoi diagram from Fig. 8(c) is a periodic framework supporting different types of stresses, (a) (Not 
F-invariant) Each of the “aligned” infinite paths supports a one-dimensional stress (equal on each edge of the path and 
illustrated here with the oblique colored diagonals). The stress values can be independently chosen on each diagonal 
path, and thus yield equilibrium stresses, such as the one depicted here, which are not F-invariant. (b) (Non-periodic 
F-invariant) A F-invariant stress assigns the same stress value on all edges in an edge orbit. Illustrated here is a stress 
where all stress orbits have the same sign; this cannot be a periodic stress, (c) A periodic stress, as the one shown 
here, must have both positive and negative stresses on edge orbits. 

We will arrive at periodic stress indirectly, via periodic liftings, defined next. In the rest of this section we show that 
the stresses corresponding to these periodic liftings must satisfy additional constraints. Then, in the next section, we 
define periodic stress in rigidity theoretic terms, and show that it coincides with this constrained stress. 



Figure 10: A periodic 3D lifting for the 
stressed framework in Figure 9(c). 

orbits of faces. Thus, there are at most n* = 


Periodic liftings. A lifting H for the planar periodic framework 
(G, F,p, n) is called periodic when F-invariant, that is, when H{p + A) = 
H{p), for all p eM? and A G A = 7r(F), where translation by periods 
has been written additively. For a face U and its translate f/ + A, this 
condition implies that: 

V[/+A = Vt/ 

Cu-Cu^x = yu’^ ( 10 ) 

Fig. 10 illustrates a periodic lifting which corresponds to the stress illus¬ 
trated in Fig. 9(c). 

Since F-invariant stresses may not always correspond to F-invariant lift¬ 
ings, additional properties are needed to characterize stresses induced by 
F-invariant liftings. We proceed now to find them. 

Stress induced by a periodic lifting. When expressed as the system 
H = (V[/,C[/)[/eF, a periodic lifting has constant coefficients Vu on F- 
card{f/r) distinct normal directions to lifted faces. 


For a closer investigation of the associated stress, we introduce the following notational conventions. A face-path 
from a face G to a face V will be indicated by G V, in particular, a face-path from G to its translate G + A will be 
indicated by G —> G A. Sums over face-paths or face-cycles are assumed to be written according to the orientation 
rule given through tetrads. With this convention, we rewrite the relations (5) and (6) obtained in Section 2 as: 


Vy = V[/+ Y, 

Cv=Cu- Y ^i\^i Pi\ ( 11 ) 

U^V 

Applying (10) we obtain the following two conditions on the stress s: 
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( 12 ) 


Y ^MiPi\=(^u-Cu+^ = Vu-X (13) 

We summarize these observations as: 

Proposition 4. If the planar non-crossing periodic framework (G,r,/7,;r) has a periodic lifting H = {Vu^Cu)u^v, 
then the associated T-invariant equilibrium stress s satisfies conditions (12) and (13) for any face U and period vector 
A G A = :7r(r). 

From constrained stress to periodic lifting. Our next goal is to show that conditions (12) on a F-invariant equi¬ 
librium stress are sufficient for determining a periodic lifting. Obviously, relations (13) will serve for identifying the 
lifting data^ = (V[/,C[/)[/eF- 

Proposition 5. Let s = {se)eeE be a F-invariant equilibrium stress for the planar non-crossing periodic framework 
(G,r,/7,;r). If for some face Uq and generators Ai,A 2 of the period lattice A = 7r(r), 5* satisfies the additional 
conditions that: 


Y Si{qi-pi)=0, 7 = 1,2 (14) 

then the lifting H = (Vjj^Cu) defined, for all A G A, by: 

V[/'A= ^ Si\qi pi\=Cu-Cu^x 

u^u+X 

is a periodic lifting inducing s, determined up to a choice of constant Co = Cuq. 

Proof Let us assume that 5* = (si)i is a F-invariant equilibrium stress which satisfies (14) for some face Go and two 
generators Ai, A 2 of the periodicity lattice A = 7r(F). We first show that our assumption implies that condition (12) is 
satisfied for all U and A G A. It is clearly satisfied for Go and all A G A since the edge vectors implicated in the sum 
are, by periodicity, the same linear combination with integer coefficients as given in A = ^lAi +^ 2 A 2 . 

Suppose next that V is adjacent to Go and consider the face-cycle Go—^y + A^Go + A->Go. Again, by 
periodicity in the framework, the edge vector implicated for Go y is the opposite of the edge vector implicated for 
y + A —> Go + A and (12) holds for V. This is enough for our claim since G and G* are connected. 

As shown in the proof of Prop. 2, for any initial choice (Vuq , Cufi) = ( Vo, Co), we obtain a lifting H = (vu^Cu) inducing 
5'. We aim at finding a periodic representative in this equivalence class. 

Let H = (Vu^Cu) be the lifting under investigation. Relations (11) will hold. Thus Vu are constant on F orbits in F. 
We show now that the function 


^u{^)= Y ^MiPi\ 

U^U^X 

is independent of the orbit representative G. Indeed, using (12) we find 

^ Si\qi-\-PiX-= 

{U >-(t/ -i-jLi)-i-A 

= Y ^MiPi\+ Y Si\qi-Piix\= Y ^i\'iiPi\ 

U^U+X U^U+X U^U+X 

Moreover, ^u{^) is linear in A and by (11): 


JFu{X) = Cu -Cjjj^x 


(15) 


(16) 
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By comparing Cy — Cu and Cy^x ~ based on (11), we obtain 


Cyj^X-Cij^X = - Y, Si\qi^X piXrX\=Cy -Cu - Y Si\qi-piX\ = 
U^V u^v 

= Cy—Cu— Y Pi)^ ' —^V—Cu — {Vy— Vu) ' 

U^V 

This gives 


^{X) = Cy —Cyj^x = Cu —Cu^X + {^V — V[/) ’X = ^{X) H- {Vy — V[/) • X 


and shows that 


^{X) — .Aif {X) = {vy — Vu)' X (17) 

This means that if we choose our initial Vq = Vuq by Vq • A = Auq{X), for all A G A, the resulting lifting will be 
periodic. □ 

In the next section we prove that this type of constrained stress is precisely the periodic stress implicated in the 
deformation theory of periodic frameworks introduced in [3]. 


4 Periodic deformations and stresses 

We have arrived at one of the most important aspects of this paper, which brings in the connection with the infinites¬ 
imal rigidity or fiexibility of a periodic framework. The application to periodic pseudo-triangulations and expansive 
mechanisms presented in the next section relies on this correspondence. 

Periodic deformations. A planar framework (G,r,/7, n) was defined by a placement of vertices p : V ^ and a 
faithful representation ;r : F —^ ^ (M?) of the periodicity group by a rank two lattice of translations A = ;r(r), with the 
necessary compatibility relation. In the framework, the edges of the graph are now seen as segments of fixed length, 
forming what is called in rigidity theory a bar-and-joint structure. According to our periodic deformation theory, 
introduced in [3] and pursued in [4, 6], a periodic bar-and-joint framework is said to be periodically flexible if there 
exists a continuous (or equivalently smooth) family of placements, parametrized by time t, : V ^ with pQ= p, 
which is not given by global rigid motions and satisfies two conditions: (a) it maintain the lengths of all the edges 
^ G E, and (b) it maintains periodicity under F, via faithful representations : F ^ (M?) which may change with t 
and give a corresponding variation of the periodicity lattice A^ = ;r^(F). 

To represent Kt we first choose two generators for the periodicity lattice F. The corresponding lattice generators X\ {t) 
and X2{t) at time t may be viewed as the columns of a non-singular 2x2 matrix denoted, for simplicity, with the 
same symbol At G GL(2). The infinitesimal deformations of the placement {pt^Tlf) are described using a complete 
set of n vertex representatives for V/F, i.e. the vertex positions are parametrized by (M^)^. The m representatives for 
edges mod F are then expressed using the vertex parameters and the periodicity matrix A. An edge representative j8 
originates in one of the chosen vertex representatives i = i{P) and ends at some other vertex representative j = j{P) 
plus some period Ac^g, where c^ is a column vector with two integer entries. The edge vectors ej^, P G E/F thus have 
the form: 


ej^ = {xjAAci^)-Xi, P G E/F (18) 

By taking the squared length of the m edge representatives, we obtain a map: 

(M^)" X GL{2) R"* (19) 

and the differential of this map at the point under consideration (i.e. the point of (R^)" x GL{2) C R^"+"^ corre¬ 
sponding to the framework (p, k)), seen as a matrix with m rows and 2n + A columns is called the rigidity matrix 
R = R{G,r,p, n) of the framework. The row corresponding to the edge described above has the form: 
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(0...0 -e'p 0...0 e'p 0...0 cje^p) (20) 

where is the transpose of the column edge vector and an obvious grouping convention is used for the columns. 

The vector space of infinitesimal periodic motions (or infinitesimal periodic deformations) of the given framework 
(G,r,/ 7 ,;r) can be described as the kernel of the rigidity matrix R and the vector space of periodic stresses can be 
described as the kernel of the transpose R^ It is understood that a stress described on the m representatives for E/F is 
extended by periodicity to all edges. 

Thus, non-trivial periodic stresses express linear dependences between the rows of the rigidity matrix R. Grouping 
these dependences over groups of columns corresponding to vertex representatives, we obtain immediately the fact 
that a periodic stress is necessarily a F-invariant equilibrium stress as defined earlier in this paper. However, there are 
two additional vector conditions imposed by the columns corresponding to the infinitesimal variation of the periods. 

This sets the stage for a comparison of the periodic stresses reviewed here and the stresses induced by periodic liftings. 
It is useful to restate now the definition. 

Definition 6. A periodic stress/ 6 >r the framework (G,F,/ 7 , n) is a stress induced from an element in the kernel of the 
transposed rigidity matrix R^ that is, a T-invariant equilibrium stress s satisfying the conditions 

L 7=1,2 (21) 

peE/r 

with integer coefficients c^^ as given in the edge description (18). 

Remarks. The fact that periodic stresses do not depend on the choices of representatives used in expressing the 
rigidity matrix follows from the fact that the image of the map (19) and the image of its differential do not depend on 
these choices. Periodic stresses give the orthogonal complement in for the image of the differential. 

There is an equivalent form for the periodic stress conditions (21), which refers directly to the periods Xp = Acj^ 
implicated in the edge descriptions (18). It can be given in terms of tensor products, namely: 

spXp®ep=0 ( 22 ) 

peE/r 

The equivalence follows from a straightforward linear algebra verification. 

For a comparison of conditions (14) and (21), we verify first the persistence of periodic stresses under relaxation of 
periodicity from F to a subgroup of finite index F C F. 

Proposition 7. Let s = be a periodic stress for the periodic framework (G, F,/?, n). Let T (ZT be a subgroup 

of finite index. Then s remains a periodic stress for the framework with relaxed periodicity (G,r,/7, 7r|p). Moreover, if 
a F-invariant equilibrium stress becomes periodic for a relaxed periodicity f C F, it is already periodic for F. 


Proof. A F-invariant equilibrium stress is obviously f-invariant. We look at the columns related to periods in the two 
rigidity matrices. 

Let p = card(r/F) denote the index of relaxation of periodicity. For the relaxed case we use the following repre¬ 
sentatives. We first choose representatives in F for the elements of the quotient group F/f. Then, a complete set of 
vertex representatives mod f will consist of the representatives chosen mod F translated by the periods correspond¬ 
ing to the representatives chosen for F/f. The subscripts for the relaxed case take the form of pairs /, 7 , with the p 
representatives 7 covering F/f. 

The proof amounts to verifying identities which can be better controlled when using a tensor product expression for 
the stress condition on the periodicity columns as in [3], page 2641 and as recalled above in (22). Altogether, with 
representatives for periodicity F given in coordinates as 

Xi, i=l,...,n 
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=Xj + i = i{^),j = j{p),Hp SA 

we have representatives for periodicity f given in coordinates as 

Xij = xi-hy 


e^y — (vj + 7) + /ijg — {xi + 7 ) — ^/3 

where, for simplicity, no further notational distinction has been made for 7 G F and its image by n. Now, we must take 
into account that we have unique expressions 


M/ 3 —7/3+'^/3 (23) 

with 7 ^ among the chosen representatives for F/f and 7 G f. Moreover, as we add representatives, we must record 
relations of the form 


7+7/3 = (24) 

with A /3 y among the chosen representatives for F/f and A /3 ^ G f . 

The periodic stress condition for F reads 

(25) 

P 


and the one for f reads 


y^sp(Xp+lpy)(S)ep=Q. (26) 

P,Y 

In order to see that these conditions are equivalent, we notice that (23) gives 


Y^splp^ep =pY^splJ.p®ep-pY^spYp®ep (27) 

P,y P P 


Similarly, from (24), we obtain 


Y^spXpj(»ep = '^sp{y+yp-?(.p^y)(S)ep. (28) 

P,Y P,Y 


Now, we take into account the identity: 


Y^sp{Y-Xpy}®ep=0 (29) 

P,Y 

which follows from the fact that when operating with an element of a group, one obtains a permutation of the elements 
in the group. Here, for any j 8 , the two lists the two lists 7 and A /3 ^ of representatives for F/f are made of the same 
elements. 

Thus, equations (27), (28) and (29) imply 

Y,^p(^P+^P,y)'^^P — pJL^pl^P'^^P ( 3 *^) 

P,Y P 

and establish the equivalence of (25) and (26). Our proposition is proven. □ 
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Comment. Tedious as it may be, this verification has the important consequence that upon relaxation of periodicity, 
the dimension of the space of periodic stresses can only go up or stay the same. 

We recall (from [3], page 2641) the relation 


G — 5 = m — 2n — 4 (31) 

connecting periodic stresses and infinitesimal deformations, where a denotes the dimension of the space of periodic 
stresses and 5 is the dimension of the space of infinitesimal periodic deformations. Subtracting the trivial infinitesimal 
deformations induced by infinitesimal isometries, we obtain: 

(T = 0-l + (m-2^) (32) 

where (j) denotes the dimension of the space of infinitesimal flexes (j) = 5 — 3. This formula will be relevant for 
evaluating behavior under relaxations, with a and (j) non-decreasing and the term (m — 2n) multiplied by the index of 
relaxation p. 

The last ingredient needed for the proof of our Main Theorem is Lemma 8 below, which can be summarized as saying 
that/6>r sufficiently relaxed periodicity f, a complete set of edge representatives can be found that are either inside a 
fundamental parallelogram Pp or cross the border with an endpoint in a neighboring parallelogram, intersecting one 
of two prescribed spanning sides. The Lemma gives a more transparent interpretation for the conditions (21) satisfied 
by periodic stresses, in terms of a sufficiently large relaxation of periodicity. 

Lemma 8. Let (G, L, p, 7t) be a planar non-crossing periodic framework. One can find generators Ai, 2.2 for the lattice 
of periods A = 7r(r) and large enough positive integers r\ and r 2 such that the relaxation of periodicity given by the 
sublattice L of index p = rir 2 of L corresponding to the two generators Xj = rjXj, (7 = 1,2) allows the following 
setup: 

(a) a choice of fundamental parallelogram Pr spanned by Ai,A2 such that its boundary avoids vertices and its 
vertices avoid the framework; 

(b) a representation of the associated fundamental parallelogram Pf for F, spanned by Xj = rjXj, (7 = 1,2) as 
tiled by p = rir 2 translated copies of the previous parallelogram; 

(c) using as representatives of vertex-orbits for L all the vertices inside the fundamental parallelogram Pr and 
then as representatives of vertex-orbits for F all the vertices inside the fundamental parallelogram Pf, a description 
of edge-orbits as used for the two rigidity matrices (for F and the relaxation F) involving, for F, all edges inside the 
parallelogram in (b) and those crossing the two spanning sides given by XjJ = 1,2; 

(d) the A periods implicated in the description of the above crossing edges are either —X\ or —X 2 or possibly their 
sum —{Xi + A2) . 



Figure 11: After a first relaxation one obtains the yellow parallelogram as fundamental domain and all edges with one 
end in it have the other end in it or in one of the eight neighbours. 

Comment. The intention, content and principle of proof of this lemma can be elucidated by referring to Figure 11. The 
aim is to show that convenient choices of lattice generators and relaxation can provide a setting where a fundamental 
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parallelogram for the relaxed lattice A is the central (yellow) zone, marked OA\{A\ H-A 2 )A 2 , with all representatives 
for edge-orbits contained in the union of the central zone with the three neighbours (in green) around the O comer. 
With this corner taken as the origin, the two generators ^jj = 1,2 of A are represented by vectors OAjJ = 1,2, we 
see that any edge in the union of the central zone and the blue zone around A i +A 2 has, by translation with — (Ai + A 2 ), 
an equivalent representative in the desired union (of yellow and green). Thus, the lemma can be proven by showing 
that convenient choices lead to a setting with no edges crossing from the central (yellow) zone to the remaining two 
(red) zones abutting corners Ai and A 2 . In fact, by translation with A 2 — Ai, it is enough to insure that no edge crosses 
from the central parallelogram to the (red) neighbor abutting at Ai. 

Proof. A generic choice of origin O will easily satisfy the avoidances desired in (a) when we represent the periodicity 
vectors as emanating from O. Assuming such a choice, rational will refer to the common rational grid determined by 
the initial periodicity lattice (and containing all its relaxations). 



Figure 12: The horizontal line h intersected by edges. The marked vector is a period along h. Since all edges have 
nearly vertical directions a sufficiently more slanted line achieves the property that no edge crosses from the yellow 
zone to the red zone. 

Since there is a finite number of edge vectors, their length is bounded and it is clear that for any initial choice of 
fundamental parallelogram Pr we can find a relaxation with all edges originating in Pf ending either in itself or one of 
the eight neighboring translated copies. Thus, what remains to be argued is how to obtain the aditional property that 
no edge crosses from the central parallelogram Py = OA\{A\ + A 2 )A 2 to the (red) neighbor abutting at Ai. 

For this purpose, we remark that our problem is not affected by a rational linear transformation and we may use 
this invariance to arrange and assume that all our edge vectors are, in direction, sufficiently close to a single direction, 
which we designate as our “vertical”. We look now at the line through O which is orthogonal to our vertical and refer 
to is as our “horizontal” line. 

Let us depict (with their periodicity along this horizontal) all the edges crossing it, as illustrated in Figure 12. Then, 
we can take a more slanted rational line with respect to the vertical as our direction for A 2 . Under our assumptions, 
this choice obtains the property that no edge crosses from the yellow region to the red region. We can find first F- 
periods along the horizontal and slanted directions. In fact, by the relative freedom we have when choosing the slanted 
direction, we may assume that we obtain generators Ay ,7 = 1,2. As already explained, a relaxation f can be found 
which satisfies the setting in Figure 11 and the additional property that all edge representatives for f can be found in 
the union of the four parallelograms around O. This completes the proof of our lemma. □ 

All the elements for the correspondence between periodic liftings and periodic stresses are now in place to prove: 

Main Theorem Let (G,r,p, n) be a planar non-crossing periodic framework. A stress induced by a periodic lifting 
is a periodic stress and conversely, any periodic stress is induced by a periodic lifting, determined up to an arbitrary 
additive constant. The correspondence relates the stress signs to the mountain/valley types of the lifted edges. 
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Proof. We use Proposition 7 and the setting described in Lemma 8 obtained after an adequate relaxation of periodicity 
f C r with generators related by Xj = rjXjJ = 1^2. We first observe that, for periodicity f the stated correspondence 
between periodic liftings and periodic stresses becomes obvious, since conditions (21) and (14) ask exactly the same 
thing: that the stress-weighted sums of edges implicated along U ^ t/ -f Ay, 7 = 1,2 be zero. The case of full periodicity 
r now follows from Proposition 7 and the corresponding fact that a F-invariant lifting which is f-periodic for some 
relaxation f C F, must be already F-periodic, as immediately seen from conditions (14). The sign relationship follows 
from Proposition 3. □ 

This concludes the proof. We turn now to applications of our Main Theorem. 


5 Periodic pointed pseudo-triangulations 

In this section we define periodic pointed pseudo-triangulations or, for short, periodic pseudo-triangulations. We 
prove properties analogous to those possessed by their finite counterparts [47, 48], with respect to expansive motions. 
In addition, we show some entirely new properties, specific to the periodic setting, related to ultrarigidity. 

A pseudo-triangle is a simple closed planar polygon with exactly three 
internal angles smaller than n. A set of vectors is pointed if there is no 
linear combination with strictly positive coefficients that sums them to 0. 
Equivalently, a set of pointed vectors lie in some half-plane determined by 
a line through their common origin. A planar non-crossing periodic frame¬ 
work (G, F, p, 7t) is a periodic pointed pseudo-triangulation when all faces 
are pseudo-triangles and the framework is pointed at every vertex. Thus, 
at every vertex, the incident edges lie on one side of some line passing 
through the vertex. As in the finite case, pointedness at every vertex is es¬ 
sential. Pseudo-triangular faces mark the ‘saturated’ stage where no more 
edge orbits can be inserted without violating non-crossing or pointedness. 
An illustration for ^ = 3 is given in Figure 13. We show that periodic 
pointed pseudo-triangulations, viewed as bar-and-joint mechanisms, sat¬ 
isfy two remarkable rigidity-theoretic properties: they have the right num¬ 
ber of edges to be fiexible mechanisms with exactly one degree of freedom 
(in the finite case [47, 48], the fiexible mechanisms were obtained after re¬ 
moving a convex hull edge), and they encounter no singularities in their 
deformation for as long as they remain pseudo-triangulations. 

Proposition 9. A periodic pseudo-triangulation has m = 2n, that is, the number of edge orbits m = card{E/T) is twice 
the number of vertex orbits n = card{\//F). 

Proof. The proof combines Euler’s formula on the torus (9) with counting the corners of pseudo-triangular faces. A 
corner is an angle smaller than n at some vertex of a given face. 

We denote by dy the degree of a vertex (valency) in the quotient multi-graph G/F, and use the classical formula 
LveG/r<^v = 2m relating the degree sum to the number of edges in a graph. Since there are dy — l comers incident to 
each vertex, we obtain a second relation: 

^ dy = n-\- 3n* , hence 2m = n-\- 3n* 
veG/r 

Combined with Euler’s formula (9) n — m-\-n* = 0 between the number of vertex, edge and face orbits, this yields 
m = 2n = 2n*. An illustration for n = 3 is given in Figure 13. □ 

Comment: Since for pseudo-triangulations we have m = 2n, formula (32) relating the dimensions (j) and cr of the 
spaces of periodic fiexes and respectively periodic stresses implies that: 

(7 = 0-1 (33) 



Figure 13: A periodic pseudo-triangulation 
with (n,m,n*) = (3,6,3). 
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Moreover, the relation remains unaffected by relaxations of periodicity of finite index. We will make use of this 
observation further down and in the next section. 

Proposition 9 shows that periodic pointed pseudo-triangulations have the right number of edges to provide smooth 
one-degree-of-freedom periodic mechanisms. The next proposition shows that this is actually the case. The argument 
amounts to showing that they have no non-trivial periodic stress. 

Proposition 10. A periodic pseudo-triangulation cannot have nontrivial periodic stresses. The local deformation 
space is therefore smooth and one-dimensional and continues to be so as long as the deformed framework remains a 
pseudo-triangulation. The same statement holds true for any relaxation of periodicity f C F of finite index. 

Proof. A nontrivial stress would give, by our Main Theorem, a periodic lifting. Such a lifting must have at least one 
vertex achieving the global maximum and another achieving the global minimum. However, with an angle exceeding 
71 on some face around each vertex, neither a maximum nor a minimum is possible. Thus, the edge length conditions 
are infinitesimally independent and the implicit function theorem gives a local deformation space which is smooth and 
of dimension one. □ 

We present now a proof of a most remarkable property of periodic pseudo-triangulations. 

Theorem 11. (Periodic pointed pseudo-triangulations have expansive Idof flexes) Let (G,r,/7,;r) be a planar 
periodic pseudo-triangulation. Then the framework has a one-parameter periodic deformation, which is expansive for 
as long as it remains a pseudo-triangulation. 

Proof outline. We compare the infinitesimal variation of the distance between two pairs of vertices (which do not 
belong to rigid subcomponents) and show that one cannot increase while the other decreases. An argument already 
used by Maxwell in [37] shows that it is enough to reason with two pairs of vertices which do not create self-crossing 
when inserted as edges (and replicated by periodicity). After one edge insertion, the framework becomes infinitesi¬ 
mally rigid and after the second edge insertion it becomes periodically stressed, with a one dimensional periodic stress 
space. Since either pointedness or non-crossing of edges is violated only by the newly inserted edges, and since a 
pointed vertex cannot be a local extremum for the height function, we infer that the associated periodic lifting must 
have its global maximum and global minimum at two of the vertices in the two added edges. Considerations of stress 
sign and of the related mountain/valley type of the lifted edges, as summarized in Prop. 3) show that the new edge 
orbits have stress factors of opposite signs. This in turn implies (by linear programming duality) that the infinitesimal 
variation of the distance between the two pairs of vertices is of the same kind: expansion or contraction. 

Proof details. We develop our arguments in a sequence of Lemmas. First, we consider what happens after one edge 
is inserted in a non-crossing manner. 

Lemma 12. The periodic framework obtained by inserting a new T-orbit of non-crossing edges in a periodic pointed 
pseudo-triangulation is minimally rigid, that is, infinitesimally rigid with 2^ -f-1 edge orbits. 

Proof. Relation (32) for the extended framework implies that the dimension of the space of periodic fiexes equals 
the dimension of the space of periodic stresses, i.e. cr = 0. We argue now that no nontrivial periodic stress exists, 
i.e. (7 = 0. Indeed, a nontrivial periodic stress would give a nontrivial periodic lifting. The latter must have a global 
maximum vertex and a global minimum vertex. Since pointed vertices cannot be either maxima or minima, the two 
extrema must occur at the endpoints of the new edge. However, a maximum vertex requires at least three non-pointed 
‘mountain’ edges and a minimum vertex requires at least three non-pointed ‘valley’ edges. The new edge would have 
to be a ‘mountain’ because of the maximum at one end and would have to be a ‘valley’ because of the minimum at the 
other end. This contradiction proves the lemma. □ 

We analyze now what happens after a non-crossing second edge orbit insertion. By inserting this new edge and its 
F-orbit in the minimally rigid framework obtained above, we must have a one-dimensional space of periodic stresses 
since the space of infinitesimal fiexes remains null. 

Lemma 13. If we add two new edge orbits to a periodic pseudo-triangulation, then we cannot have the same sign for 
the stress factors on the two new edge orbits inserted in the pseudo-triangulation. 
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Proof. We argue first for the case of no common vertex mod F for the two new edges. Then, as shown above in 
Lemma 12, we cannot have the maximum and the minimum at the two ends of the same new edge. Thus the edge 
reaching to the maximum is a ‘mountain’ and the other one, reaching to the minimum, is a ‘valley’. 

In case the two new edges share a vertex, we note that at least one extremum must be at some unshared vertex. But 
then the other extremum cannot be at the shared vertex and must be at the other unshared vertex. Thus, as above, the 
two edges must be a ‘mountain’ and a ‘valley’. □ 

With opposite signs confirmed for the stress factors on the two new edges, we look now at the infinitesimal displace¬ 
ment induced by the infinitesimal deformation of the pseudo-triangulation on these same edges. 

For the framework with two new edge orbits, the non-zero stress vector must be orthogonal on the image of the 
corresponding rigidity matrix. In particular, it must be orthogonal on the image vector obtained by evaluating the 
rigidity matrix on the infinitesimal displacement induced by the pseudo-triangulation. But all entries corresponding to 
‘old’ edge orbits, that is, edge orbits in the pseudo-triangulation, will be zero. With opposite stress signs for the two 
new edges, orthogonality requires the same sign for the two non-zero entries. This means: simultaneous infinitesimal 
expansion or simultaneous infinitesimal contraction. 

By now, the following result has been established. 

Proposition 14. Let (G,r,/7, 7t) be a planar pseudo-triangulation. Let two pairs of vertices be such that both corre¬ 
sponding distances vary infinitesimally under the infinitesimal deformation of the pseudo-triangulation. Assume that 
inserting the corresponding two edges (and their F orbits) does not produce self-crossing. Then, the two edges vary 
infinitesimally in the same way: both expand or both contract. 

As a final step to complete the proof of The¬ 
orem 11, we want to remove the assumption 
that the two new edges maintain the non¬ 
crossing nature of the periodic framework, 
and thus to prove the full expansive character 
of the pseudo-triangulation mechanism. In 
fact, because we are comparing the stresses 
and infinitesimal fiexes of pairs of edges, it 
suffices to assume that the first one is not 
crossing the rest, hence we only need to re¬ 
move the non-crossing assumption for the 
second inserted edge. 

Figure 14 illustrates the (descending) induc¬ 
tive step whereby stress considerations in a 
situation where an edge intersects several other edges can be reduced to an equivalent situation without self-crossing. 
In [37] the procedure is called Bow’s method. 

We briefiy review the nature of the argument by explaining the stress correspondence for the step shown in Figure 14. 
A new vertex has been inserted at the last crossing of the marked edge vectors. With obvious notations (by which an 
edge a is split into two edges a\ and a2, and similarly for the edge j3), we have, for the corresponding edge vectors: 

^/3 = ^/ 3 i + ^^2 

Geometric stress conditions at the new vertex require 

^ai^ai = ^a2^a2 

and the correspondence of stresses (valid for periodic stresses as well) is expressed through the formulae: 

— ^a2^oc2 
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with all other stress factors (on edges different from a and j8) remaining the same. 

Relations (34) show that in this isomorphic correspondence of stresses, the sign of the stress factors along the frag¬ 
mented edges is the same as in the initial framework, hence the argument given in Lemma 13 carries over. Thus, 
whether crossing or not crossing other edges, the insertion of a second new orbit of edges in a periodic pseudo¬ 
triangulation produces a stressed framework with opposite signs for the stress factors on the two added edge orbits. As 
shown, this implies the same type of infinitesimal variation of the two corresponding vertex distances when deforming 
the periodic pseudo-triangulation. The limit case, when the inserted edge passes through one or several vertices of the 
framework is treated similarly. In fact, the step for crossing through a vertex is simpler and only requires the splitting 
of the inserted edge into two edges. This concludes the proof of Theorem 11. □ 

Comment. We retain from this section the conspicuous property of periodic pseudo-triangulations of maintaining the 
same deformation space under arbitrary relaxations of periodicity of finite index. We make use of this property in the 
final section, where we uncover the ultrarigidity of a family of frameworks obtained from pseudo-triangulations. 

We conclude the paper with two applications motivated by questions arising in materials science. 

6 Expansive and auxetic paths 

Periodic pseudo-triangulations determine a large, yet distinctive class of planar periodic frameworks. We have pre¬ 
sented above some of their remarkable properties: their local deformation space is smooth and one-dimensional and 
this deformation path is expansive as long as the framework remains a pseudo-triangulation, that is, for the proper 
sense of variation of the parameter, all distances between pairs of vertices increase or stay the same. Moreover, a 
relaxation of periodicity to a subgroup of finite index does not change these characteristics. 

In this section we emphasize, from a purely geometric perspective, the relevance of periodic pseudo-triangulations and 
their expansive paths for what has been termed auxetic behavior in materials science. 

In materials science, the term auxetic, suggestive of increase or growth, refers to solids with a negative Poisson’s ratio 
[32, 21]. Quoting from [33]: ‘A negative Poisson’s ratio in a solid defines the counter-intuitive lateral widening upon 
application of a longitudinal tensile strain.” For anisotropic solids, such as single crystals, “the variation of elastic 
moduli with direction is also relevant” (ibid, p.6445). Thus, auxetic behavior is primarily defined in terms of physical 
characteristics of the material under consideration. However, in many instances with pronounced geometric structural 
underpinnings, geometrical explanations have been proposed [51, 24, 25, 39]. 

In our context, which is that of periodic frameworks and their deformations, we must rely on a comprehensive but 
strictly geometrical concept of auxetic path. This concept is introduced in our companion paper [8] and will be briefiy 
reviewed here. 

The geometric approach to auxetics presented in [8] is formulated in arbitrary dimension d and addresses one- 
parameter deformations of a given periodic framework in W^. Two aspects may be emphasized from the start: (i) 
that when present, the auxetic property refers to the deformation path under consideration and not the framework it¬ 
self, which oftentimes allows other deformation paths which need not be auxetic and (ii) the auxetic character is an 
expression of a particular type of variation of the periodicity lattice along the deformation path. 

Mathematically, we rely on the notion of contraction operator in a Hilbert space, which we recall below. For our 
purposes, with the standard norm |v| = will sufffice. Let T : ^ be a linear operator. Then, the 

operator norm, or simply the norm of T is 


\\T\\=sup\^\<i\Tx\=sup\^\=i\Tx\ 

T is called a contraction operator, or simply a contraction, when 11T11 < 1. 

Comment. From \Tx\ < | |r| | • |v| it follows that contraction operators are characterized by the property of taking the 
unit ball to a subset of itself. For a strict contraction one requires 11T11 < 1. 

Let us assume now that we have a one-parameter deformation (G,r,/7T, TT^), T G {—£,£) of a periodic framework in 
W^. The corresponding periodicity lattices = 7rT(r), offer by themselves a way to compare any two sequential 
moments Ti < T2 by looking at the unique linear operator 7^2 Ti defined by 
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( 35 ) 


Definition 15. A differentiable one-parameter deformation G (—£,£) of a periodic framework in 

is called an auxetic path, or simply auxetic, when for any Ti < T2, the linear operator defined by (35) is a 
contraction. 

In [8] we prove the equivalence of this geometric criterion for auxetic paths with the following characterization involv¬ 
ing the evolution of the Gram matrix of a generating basis for the period lattice of the framework. With conventions 
already used in previous sections, after choosing an independent set of generators for F, the image 7rT(r) is completely 
described via the J x J matrix with column vectors given by the images of the generators under The associated 
Gram matrix will be 


ffl(T) = VAt. 

Proposition 16. A deformation path (G,r,;7^, ;rT), t G (—£,£) is auxetic if and only if the curve of Gram matrices 
(o{t) defined above has all its tangents in the cone of positive semidefinite symmetric d xd matrices. 


Remarks. This infinitesimal characterization of auxetic paths is easily seen to be independent of the choice of gener¬ 
ators used for obtaining the Gram matrices. There is an obvious analogy here with causal lines in special relativity. 
Causal paths must have all their tangents in the light-cone of Minkowski space-time. Likewise, auxetic paths must 
have all their tangents ^(t) in the cone Q.{d) of positive semidefinite symmetric matrices. For the geometry and 
linear symmetries of this cone see [26, 1]. 


This criterion is particularly convenient when frameworks and deformations are described as in [6] using parameters 
in X Q,{d). The coordinates in describe the position of (^ — 1) chosen representatives for vertex 

orbits relative to the periodicity basis A which can be retrieved, up to orthogonal transformations from its Gram matrix 
CO = A^A G The orbit of a first vertex representative is identified with the periodicity lattice. 


■a 

Figure 15: (Left) The Kagome framework. (Right) 
Parametrizing the deformation of the Kagome framework. 
Triangle OAB is fixed, and triangle OCD rotates with an 
angle 6 from the standard position. The two generators of 
the periodicity lattice are marked as arrows. 


We now address the comparison of expansive and aux¬ 
etic paths. 

Theorem 17. Let (G,r,;7^,;r^),^ G (—£,£) be a one- 
parameter deformation of a periodic framework in 
If the path is expansive, that is, if the distance between 
any pair of vertices increases or stays the same for in¬ 
creasing t, then the path is also auxetic. However, aux¬ 
etic paths need not be expansive. 

Proof. As emphasized earlier, the auxetic property de¬ 
pends only on the curve (o{t) and it will be enough to 
use the expansive property on one orbit of vertices. We 
have to verify that the operator 7 t2Ti which takes the pe¬ 
riod lattice basis A^2 period lattice basis A^^ is a 

contraction for T2 > Ti 


Let us observe that in the unit ball of the vectors with 
rational coordinates relative to the basis A^^ ^ dense 



subset. Since some integer multiple of such a point is a period at moment T2, and this period, as a distance between two 
vertices in a vertex orbit, can only decrease or preserve its norm when mapped by Ttiti to the corresponding period 
at moment Ti, we see that a dense subset of points in the unit ball must be mapped to the unit ball. This is enough to 
conclude that 11T11 < 1. □ 


The fact that small auxetic deformations need not be expansive is to be expected from the fact that the relative motion 
of different vertex orbits is not sufficiently constrained by the auxetic property. A simple example is offered by the 
Kagome framework in dimension two. The Kagome framework is a familiar planar example and has been explored 
from various points of view [27, 28, 31, 19, 49]. The auxetic character of its one-parameter deformation is frequently 
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mentioned [24, 38, 39]. The brief review here is meant to offer a simple illustration of our geometric criterion for 
auxetic paths and to distinguish the expansive portions of this deformation. 

Example 2. (Deforming the Kagome framework) The basic elements for describing the framework are shown in 
Figure 15. The parametrization for the deformation is described in Figure 15. Triangles OAB and OCD are assumed 
congruent and equilateral. With origin at O, coordinates may be chosen so that A = (—1,0) and B = —1/2(1, \/3). 
The resulting Gram matrix for the marked periods and a rotation of triangle OCD with angle 6 from the standard 
position is: 


fi)(e) = (i+cose) ( J 2 ) (36) 

with: 

^(0) =-sine (5 J) <37) 

The image of the deformation path in the positive semidefinite cone fl(2) for 0 G (—TT, n) is the segment from the null 
matrix to co{0), covered twice. The standard position 0 = 0 is a maximum for the area of a fundamental parallelogram. 
Proposition 16 shows immediately that the symmetric paths obtained for 6 running from 7i to 0, respectively — ;r to 0, 
are auxetic. 

The range of 6 corresponding to periodic pseudo-triangulations is the union (—2;r/3, —7rl3)[J {7113^27113). Expansive 
behavior cannot occur beyond this range, since the distance variation between the two pairs of vertices (A,Z)) and 
(sS,C) has opposite character when 6 is in the complement: when one segment increases, the other one decreases. 

The framework shown in Fig. 5 is one of the emblematic illustrations of auxetic behavior in dimension two. It is called 
a reentrant honeycomb in [42]. It is often used to illustrate auxetic behavior: a vertical stretch involves a necessary 
horizontal expansion. 

Example 3. (The ‘‘reentrant honeycomb”) The framework in Fig. 5 has two degrees of freedom and not all defor¬ 
mations paths are auxetic. The expansive possibilities can be explained in terms of the two possible refinements to 
periodic pseudo-triangulations shown in Fig. 16. 

The framework with all edges of the same length 
may be obtained by deforming the 2-diamond 
framework (the regular hexagonal structure) [9]. An 
auxetic path would result from a ‘vertical stretch’ 
which visibly entails a horizontal expansion as well. 

However, an understanding of all expansive in¬ 
finitesimal deformations within all auxetic possibilities requires some elaboration [10]. Here, we limit our discussion 
to the following remarks on the basic role of periodic pseudo-triangulations. 

The local deformation space of the framework is smooth 
and two-dimensional. Indeed, by pointedness at every ver¬ 
tex, there can be no periodic stress. Since n = 2 and m = 3, 
there are two degrees of freedom. Figure 16 marks the two 
angles which may be used as parameters. By maintaining 
the central symmetry of the hexagon, the two depicted pairs 
of periods remain equal vectors. The framework can be 
refined in two ways to a periodic pseudo-triangulation by 
insertion of an additional orbit of edges. It can be shown 
that the expansive infinitesimal deformations determined 
by these two pseudo-triangulations give the extremal rays 
of the infinitesimal expansive cone of the framework. For 
the larger auxetic cone, we refer to [8]. 




Figure 17: For the ‘reentrant’ structure of hexagons in Fig¬ 
ure 5, which has two degrees-of-freedom, the figure on the right 
marks two angles which may be used as parameters. By main¬ 
taining the central symmetry of the hexagon, the two depicted 
pairs of periods remain equal vectors. 



Figure 16: Two possible refinements to periodic pseuUo- 

triangulations of the ‘reentrant’ structure of hexagons in Figure 5. 
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7 Ultrarigidity 

In this final application we discuss the relevance for ultrarigidity of the other remarkable feature of periodic pseudo¬ 
triangulations, namely the invariance of their local deformation space under arbitrary relaxations of the periodicity 
group to subgroups of finite index f C F. 

The deformation theory of periodic frameworks, as founded and developed in [3, 4, 6], is formulated for periodic 
graphs understood as pairs (G,r): the periodicity group F is an essential structural part of the data. The study of 
phenomena involving changes in periodicity or asymptotic behavior under successive relaxations of periodicity may be 
intricate. We have recently proposed [5] the notion of ultrarigidity for expressing the property of a periodic framework 
(G.F,p, 7t) of being infinitesimally rigid and remaining so under any relaxation of periodicity to f C F of finite index. 

Ultrarigidity is easily recognized in frameworks made of some rigid finite 
parts which are rigidly connected between themselves. However, this is not 
the general case, as the following constructions, based on planar periodic 
pseudo-triangulations, will demonstrate. 

Let (G, F, p, ;r) be a periodic pseudo-triangulation. Then, the local deforma¬ 
tion space is a smooth curve. The infinitesimal deformation corresponding 
to this one degree of freedom mechanism must induce on some pairs of ver¬ 
tices a non-trivial infinitesimal variation of length. Then, as argued above 
in Lemma 12, by selecting such a pair of vertices and by inserting the cor¬ 
responding orbit of edges, we obtain a minimally rigid framework, that is, 
an infinitesimally rigid framework with m = 2^ + 1. 

In fact, the resulting framework is ultrarigid. Indeed, for any relaxation 
of periodicity to a subgroup f G F of finite index, the older framework 
with relaxed periodicity remains a pseudo-triangulation and has the same 
local deformation space. Thus, the distance between the selected pair of 
vertices varies infinitesimally and the same argument applies, showing that 
insertion of the corresponding f edge-orbit already yields an infinitesimally 
rigid framework. This concludes the proof of the following result. 

Proposition 18. Let (G,F,/7,;r) be a periodic pseudo-triangulation in 
which we consider a pair of vertices with a non-trivial infinitesimal varia¬ 
tion in distance under the infinitesimal deformation of the framework. Then 
the insertion of the corresponding F orbit of edges results in an ultrarigid 
framework. 

Thus, periodic pseudo-triangulations and insertion choices provide endless 
examples of ultrarigid frameworks. 

Example 4. (An ultrarigid periodic framework) The framework illustrated in Figure 18 is ultrarigid. The colors 
of the vertex orbits indicate the periodicity lattice of the pseudo-triangulation used in this construction. The pseudo¬ 
triangulation itself is a deformed 2-diamond framework with an additional edge orbit [9]. The edge orbit which turns 
the periodic pseudo-triangulation into an ultrarigid framework creates the rigid quadrilaterals shown in the picture. 

It is worth remarking that a finite fragment of the framework may be covered by a fragment resembling the one depicted 
in the figure. Such fragments, as finite linkages, are flexible and can accommodate small variations of the segments 
A\A 2 , A 2 ,A 3 etc. However, larger assemblies of these stacked rows of rigid quadrilaterals will have smaller leeway 
for variation of deformation parameters. In the limit, as an infinite periodic framework, the structure is rigid. 

In conclusion, we anticipate that our periodic version of Maxwell’s Theorem and the expansive nature of periodic 
pseudo-triangulations will find, like their finite counterparts, further applications in discrete and computational geom¬ 
etry. In the larger scientific context, applications are expected in new materials and mechanism design. 

An extended abstract of this work has appeared in [7]. 



Figure 18: An ultrarigid framework 
obtained from a periodic pseudo¬ 
triangulation with an additional edge 
orbit. When ignoring periodicity, finite 
fragments are fiexible. 
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